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Abstract 

We calculated the color-octet contribution to the double spin asymmetry of 
J/ ip hadroproduction with nonzero transverse momenta at fixed target energies 
^/s ~ 40 GeV. It is shown that the color-octet contribution is dominant in 
the asymmetries. The expected asymmetries and statistical errors in a future 
option of HERA with longitudinally polarized protons at HERA-A^ should allow 
one to distinguish between different parametrizations for the polarized gluon 
distribution in the proton. 
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1 Introduction 



The presently most accurate way to measure the polarized gluon distribution function 
in the nucleon is to study those processes which can be calculated in the framework of 
perturbative QCD (PQCD), i.e. for which the involved production cross sections and 
subprocess asymmetries can be predicted. One of the cleanest ways to probe QCD 
is to investigate heavy quarkonia production processes. Heavy quark pair production 
processes can be controlled perturbatively due to large mass of constituents. On the 
other hand, heavy quark systems are mainly produced in gluon fusion processes and 
therefore, asymmetries are expected to be sensitive to the polarized gluon distribution 
in the nucleon. Investigation of heavy quarkonia production processes in polarized 
experiments would also yield additional information about the quark-antiquark pair 
hadronization phase. 

The two-spin asymmetry in J/ip production has been studied in the framework 
of the so called color singlet model (GSM) by Morii and collaborators 0. But as 
was shown in the last years, the color singlet model does not describe satisfactorily 
the heavy quarkonium hadropro duct ion at Tevatron and also at fixed target energies. 
While the CSM gives a reasonable description of the J/ip production cross section 
distribution shapes over px or xp at fixed target energies, it completely fails in the 
explanation of the integrated cross section (a K factor 7-10 is needed to explain 
experimental data) . The anomalously large cross section 0] of the J/ ip production 
at large transverse momenta at the Tevatron revealed another negative feature of the 
CSM. Within the framework of the CSM it is impossible to explain the anomalously 
large ip @ and direct J/ip production in the CDF experiment at the Tevatron. 

The CSM is a nonrelativistic model where the relative velocity between the heavy 
constituents in the bound state is neglected. But discrepancies between experimental 
data and the CSM predictions hint that 0{v) corrections as well as other mechanisms 
of quarkonium production, which do not appear in the leading order in v, should 
be considered. Expansion of quarkonium cross sections and decay widths in powers 
of relative velocity v of heavy quarks in a bound state has recently been realized in 
terms of Nonrelativistic QCD (NRQCD) This formalism implies not only color- 
singlet processes but the new color-octet mechanism, when a quark-antiquark pair 
is produced on small time scales in color octet states and evolves into a hadron by 
emission of soft gluons. The color octet mechanism takes into account the complete 
structure of the quarkonium Fock space while in the CSM only the dominant Fock 
state is considered, which consists of a color singlet quark antiquark pair in the definite 
angular-momentum state (higher order Fock states are suppressed by powers of v). 
According to the factorization approach based on the NRQCD, the production cross 
section for a quarkonium state H in the process 

A + B ^ H + X (1) 
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can be written as 



^ X] / dxidx2fi/A{xi)fj/B{x2)criij H) (2) 
a{z3^H)=Y.C'^[nmO"[nm 

n 

where jijA is the distribution function of the parton i in the hadron A. The sub- 
process cross section is separated into two parts: short distance (C*-^ [n]) coefficients 
and long distance matrix elements (0|O^[n]|0). The C'-' fn] is the production cross 
section of a heavy quark- ant iquark pair in the i and j parton fusion. It should be 
calculated in the framework of pQCD. The [n] state can be either a color singlet or a 
color octet state. The (0|O^[n]|0) describes the evolution of a quark- ant iquark pair 
into a hadronic state. These matrix elements cannot be computed perturbatively. 
But the relative importance of long distance matrix elements in powers of velocity v 
can be estimated by using the NRQCD velocity scaling rules 0. 

Shapes of the distribution of short distance color-octet matrix elements indi- 
cate that the new mechanism can explain the Tevatron data of direct Jjil) and ■?/'' 
production at large j9t @- But unlike color-singlet matrix elements connected with 
the subsequent hadronic nonrelativistic wave functions at the origin, color octet long 
distance matrix elements are unknown and should be extracted from experimental 
data. The color octet contribution to the J/ t/' photoproduction has been analyzed in 



the papers [jTO|, [TT|. Recently, the J/V' hadroproduction at fixed target energies has 



been studied by including the color-octet mechanism |T^, [1^. Large discrepancies 
between experimental data and the CSM predictions for the total cross section of 
the Jj^l) hadroproduction were explained. The color octet contribution is dominant 
in the J j^^ hadroproduction at energies -^/i ^ 30 — 60 GeV. The analyses carried 



out in these papers |T0|, |TT], [T^, |13[ demonstrate that fitting the photoproduction and 
hadroproduction data at low energies requires smaller values than those extracted 
from prompt J/^^ production at CDF at the Tevatron 0. Possible reasons for such 



discrepancies have recently been analyzed in the papers |T^, |T1|. The extraction of 
color-octet long distance matrix elements from the quarkonium production properties 
in polarized experiments would be an additional test of universality. 

In the present letter we consider the color octet contribution to the double spin 
asymmetry of J j^l) hadroproduction. Unlike the previous calculations |jl6| we consider 
(cc) color octet and color singlet pair production in 2 ^ 2 subprocesses to obtain 
asymmetries at nonzero transverse momenta (p^ > 1-5 GeV). Such values of can 
not be caused by internal motion of partons in the nucleon and hence transverse 
momentum distributions of production cross section and asymmetries are calculable 
perturbatively. The double spin asymmetries in parton collisions are presented in 
section 2. Since heavy quarkonium is mainly produced in the gluon-gluon fusion sub- 
processes, the J/V^ production asymmetry should be sensitive to the polarized gluon 
distribution function in the proton. We have calculated the expected asymmetry of 



3 



J/ ip production at HERA- A^, one of the future options of HERA ; an experiment 
utilizing an internal polarized nucleon target in the polarized HERA beam with en- 
ergy 820 GeV would yield -y/i ^ 39 GeV. For comparison, we also considered the 
expected asymmetry of J/ip production in similar spin physics experiments at much 
higher energies at the RHIC collider [jlj 



2 Double Spin Asymmetries at Subprocess Level 

Let us discuss the two-spin asymmetry All for the inclusive J/ip production which 
is defined as 

j/^ ^ da{p+p+ ^ J/tjj)- da{p+p^ ^ J li^) ^ EdAa/d^p 

^^^^ da{p+p+ ^ J/V') + da{p+p^ J/^) Eda/d^p ' ^ ' 

where denotes helicity projection sign on the proton momentum direction. 

The cross section of J/ip production can be written as 

cxj/v- = cT{J/^)dir + J2 Brixcj J/^X)a^^j + Br{^lj' Jl^X)a^, (4) 

J=0,l,2 

where Br{^[cc) — > J/ipX) denotes the branching ratio of the corresponding (cc) state 
into J/ip. The production of each state of quarkonium is contributed by both color 
octet and color singlet states, as in the case of direct J/ip production 

a{Jl^U = o-f + a^ji^ = a(J/^)-"^' + Y.^iQQ['^^'L'j]moi/^C'^'Lj)\0) (5) 

where the sum is over the states ^Pqi 2, and ^Sf. We consider only the dominant 
sets of color octet states by the NRQCD velocity expansion for the direct S and P 
state charmonium production. 

In a recent paper was considered asymmetry of J/ip hadroproduction in the 
color octet model exploiting only the lowest order subprocesses (2 — (cc)) over the 
QCD coupling constant. Subprocesses 2 — 1 contribute only to the production of a 
quarkonium state at zero transverse momentum with respect to the beam axis. Trans- 
verse momentum distributions of (cc) states and, consequently, the J/ip meson are 
not calculable in the pt < ^qcd region. To deal with the experimentally observable 
quantities (taking into account the kinematic restriction on the angle with respect to 
the beam axis at HERA-A^ and at collider experiments), we consider the subprocess 
2^2 which gives the leading contribution to the quarkonium production with pt 
greater than 1.5 GeV. Such large transverse momenta cannot be caused by internal 
motion of partons in the nucleon and, respectively, subprocesses 2^1 should not 
contribute to the production of quarkonia with px > 1.5 GeV. For calculating the 
expected asymmetries we consider only the following subprocesses: 

9 + 9 (cc) + g 

g + q -> (cc) + q (6) 
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Figure 1. Partonic level double spin asymmetries for ^5o and ^Si color octet and ^Pq color 
octet and singlet states versus tj and xt- 

for the color octet and singlet states of a (cc)-pair. The quark-antiquark collision 
subprocesses are not taken into account, because of small sea quark polarization in 
the proton. 

Using the symbolic manipulation program FORM , we have calculated the (cc)- 



states production cross sections for different helicity states of colliding partons. For 
the color octet states total cross sections have been calculated by Cho and Leibovich 
p. The total cross sections for color singlet states are presented in The 
latter results served for checking our calculations. 

Figs.l and 2 present the values of the subprocess level asymmetries for the color 
octet and color singlet states production of a (cc)-pair depending on the two dimen- 
tionless quantities: r] = Aml/s^ and Xt = Vr/Vmax^ where Pmax is the maximum 
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Figure 2. Partonic level double spin asymmetries for ^Pi and '^P2 color octet and singlet 
states versus r] and xt- 

momentum of the produced state in the subprocess. In Figs.l and 2 are shown only 
the gluon fusion subprocess asymmetries because they give the main contribution to 
the hadronic level asymmetry. The cases of ^5*0 and ^Si singlet states are omitted 
in figs.l and 2. The ^Sq state does not contribute to J/ip production and analytical 
expressions for the corresponding cross sections for the ^Si state are given in 0. 

As can be seen from figs.l and 2, the asymmetries are very similar for color octet 
and color singlet states with the same spin-orbital quantum numbers. It is worth 
mentioning that in the limit 4m^/s 1 (i.e. at the threshold of heavy quark pair 
production) qll for the scalar and tensor states tends to 1 and -1, respectively. The 
limit 4m^ — > s means that the emitted gluon in the 2^2 subprocess {gg — > {cc)g) 
becomes soft and the helicity properties of the amplitude should be the same as 
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those of the amphtude of the 2 — >• 1 process {gg (cc) 0. It is easy to show 
that the asymmetries for the process 2—^1 for scalar and tensor states are 1 and 
-1, respectively. The existence of such limits serves as an additional test for our 
analytical calculations of cross sections Aa. for scalar and tensor states are 1 and -1, 
respectively. The existence of such limits serves as an additional test for our analytical 
calculations of cross sections Aa. 



3 Matrix Elements 

For the calculation of the hadronic level asymmetries of J/ production we used the 
long distance matrix elements fitted from various experimental data. All color singlet 
matrix elements are related to the radial quarkonium wave functions at the origin and 
their derivatives. As in paper p| for this purpose we used the Buchmiiller-Tye wave 



functions values at the origin tabulated in ref.[jT3[. The color octet matrix elements 
were fitted from the J/ip and higher charmonium state production data in various 
experiments. Unfortunately, there are some discrepancies between the values of color 
octet matrix elements extracted from different experiments. 

The number of color octet long distance matrix elements should be reduced by 
using the NRQCD spin symmetry relations: 

{0\O^CPjm = (2J+l)(0|OfePo)|0), (7) 
{0\Or'CSi)\0) = {2J+l){0\OreS,)\0)- (8) 

These relations are accurate up to v"^. 

After using these relations we have only four independent matrix elements 
(Os^^i^Si)), (Cf {Oi^'^i^'Po)) and (Oi^^CSo)), which give the main contri- 

butions to the J/tp hadroproduction cross section. For the first two parameters we 
used the values extracted from Tevatron data by Cho and Leibovich 0. The other 
parameters, connected with the ip' production, were also taken from [P]. As con- 
cerns the parameters {Og^^{^So)) and {Og^^ {^Pq)) , it is possible to extract only their 
combination from direct J/ip production data at Tevatron [Q: 

{OlOi^^CSom + A,{0\Oi/^ePom = Q.Q ■ lO-'CeV^ (9) 

Values for the other combinations are extracted from the J/ijj photoproduction and 



fixed target hadroproduction data [|T0], |ri|, [T^, [T^ 

{0\Oi^\'So)\0) + l-{0\oi/^{^P,)\0) = 2 ■ 10-'GeV'[13], 



{0\Oi^^CSom + ^(O|08^/^(='Po)|O) = 3 ■ HI- (10) 



■^The same reason makes the asymmetry less sensitive to the radiative corrections, than spin- 
averaged cross section, as K-factor is dominated by the contributions of virtual, soft and coUinear 
gluons 
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If one assumes that (0|C8^'^(i5o)|0) = (<d\Oi'^{^Po)\Q) /ml the photoproduction and 
fixed target hadroproduction values ('low energy' values) are an order smaller than 
the Tevatron value ('high energy' value). The possible sources for such a discrepancy 
are discussed in |jl3|, |I4| . In paper [0] it is mentioned that the mass of the produced 
hadronic final state (or intermediate color octet state) must be higher than that of 
the J /if) meson - Mj/^ = 2mc, because the intermediate octet state emits gluons with 
energy 2mcf ^ before transition into a color singlet state. But a charmonium is not a 
truly nonrelativistic system, the average relative velocity of constituents in J / ip is not 
very small — ~ 0.23 — 0.3 and emitted gluons in the cc pair hadronization phase 
have the energy 0.7^ GeV. Therefore the mass of the hadronic final state is about 
4 GeV. At fixed target energies, increasing the mass of the produced hadronic state 
leads to a reduced partonic luminosity because of a steeply falling gluon distribution 
function. Consequently, the production cross section of the (cc) color octet state 
should be smaller. So, the 'true' color octet long distance matrix elements must be 
larger, than those extracted by using Mj/^ as mass of the intermediate octet state 

Another possible source was mentioned in the paper and should be an un- 



certainty connected with the choice of different parametrizations for the gluon dis- 
tribution function in fitting the direct J/ip production data at CDF 0. Cho and 
Leibovich used in their calculations the MRSDO parametrization for parton distribu- 
tion functions. At small values of pt — 5 GeV, using a more reliable parametrization 
for small partonic x (GRV LO or GRV H0[^) leads to 1.5 -r- 1.6 times higher cross 
sections for the ^Pq and ^5*0 color octet state production than those obtained by us- 
ing MRSDO parametrization 0. Hence the fitted value of combination (|^) should be 
approximately 1.5 -^ 1.6 times smaller: 



Oi^^CSom + —{0\Oi^'^CPom = 4 - 4.4 ■ 10-^GeV'[U]. (11) 



m 



At large — 10 ^ 15 GeV all parametrizations give practically the same magnitude 
for the color octet states cross sections. Hence if we use the new combination ([Tip 
one needs a larger value of the parameter {0\Oq^'^{^Si)\0) to explain the experimental 
data for the J/ip cross section at pt — 10 ^ 15 GeV: 

(0108/^^(351)10) ~ 10 • 10-''GeV\ (12) 

Practically the same values for these parameters were obtained by Beneke and Kramer 
in a recent paper by fitting the CDF data of direct J/ip production and using the 
GRV LO (1994) parametrization 



(0 1 08/^^(351)10) = 1.06 ± O.Utlil ■ (13) 
{0\Oi^'^CSom + ^(0|O8^/^('n)|0) = 3.9 ± 1.15i}:^? ■ lO'^GeV'. (14) 
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For calculations of the liadronic level asymmetries we use three different sets for 
the three long-distance color octet matrix elements 

A) - {0\Oi^^CSi)\0) = 6.6 ■ 10-=^GeV-3[9], 
{0\Oi/\'So)\0) + A(o|Os^/'^(3po)|o) = 6.6 • lO-'GeV' 

B) - (0108/^^(^51)10) = 6.6 ■ lO'^GeVS^ 
{OlOi^^'Som + A(0|08^/'^(3Po)|0) = 3 • lO-^Gey^ 0; 

C) - {0\Oi^'^CSi)\0) = 10 ■ 10-=^GeV-3H, 
{0\Oi^\'So)\0) + A(0|08^/'^(3Po)|0) = 4 ■ lO-'GeV' (0). 

As has been mentioned above, values for the other color octet parameters were taken 
from IP . For the calculations of the expected asymmetries we assume that the leading 
term in the combination (|Tl|) is the parameter {Og^^{^So)) (or (Og'^'^(^Po)) = 0). Only 
wiyh such a radical choice the values of 'high' and 'low' energy parameters should be 
consistent to each other. 



4 Results and Discussion 

Fig. 3a presents the expected asymmetries for the three sets of color octet parameters 
in the case when the parameter {Og^^ {^Pq)) tends to zero {{Oi^'^i^So)) is the leading 
term in the combinations @ , (|10D , (|ll]) ) . Fig. 3b represents the other radical choice, 
when the second parameter {O^^^ Sq)) is zero. For the calculations of hadronic 
level asymmetries we used the GRV LO parametrization for unpolarized distribution 
functions and a parametrization proposed by Gehrmann and Stirling for the 
polarized parton distribution function (set A) For the mass of the charm quark 



the value rric = 1.48 GeV was taken. As we can see from Fig.3, the magnitude of the 
asymmetry is more sensitive to the choice of the leading term in the combinations than 
to the choice of the particular set of the octet matrix elements. Recently the color 
octet matrix elements {O^^^ {^Pq)) and {Ol^^ (}Sq)) have been extracted separately 
from the J/^ electroproduction data p^ . 

{Q\Oi''^{^SQ)\Q) = A-lQ-^GeV\ 
mj. 

These values confirm our assumption that the leading term in the combinations is the 
parameter {O'l^^ {}Sq)) . In fig. 3a, the dotted-dashed curve corresponds to the asym- 
metry calculated by using the values of the long-distance parameters (pIS]). For the 
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4 4.5 

PT(GeV) 



Figure 3. The double spin asymmetries for different sets of color-octet long distance param- 
eters. Solid line corresponds to set C, dashed line - set A, dotted line - set B; (a) for the 



case when {Q-^^"^ Pq)) 



0, (b)for - (Og^'^(^5o)) = 0. Dot-dashed line in Fig.a corresponds 



to asymmetry calculated by using values of parameters from 



parameter {Og^^ {^Si)) a value from set C was used. As one can see from figs.3(a,b), 
the expected asymmetries for all sets of parameters are practically the same. In 
the following calculations we shall use the set C for the values of the parameters, 
assummg {0\Oi^^{^Po)\0) = 0. Fig.4a shows the expected double spin asymmetries 
at HERA-A^ energies as functions of transverse momentum J/ip in the c.m.s. The 
solid curve corresponds to the Gehrman and Stirling parametrization for polarized 
distribution functions, set A; the dashed curve, to set C [^. In fig. 4a we also display 
the expected statistical errors SAll at HERA-A^ which can be estimated from 



SAll = 0.17/ J aipb) 



(16) 



This relation has been determined by assuming an integrated luminosity of 240 pb~^ 
and beam and target polarizations Pb = 0.6, Pt = 0.8 [jl5|. The error bars are ob- 
tained by using integrated cross sections over bins Apx = 0.5 GeV (for the first three 
points) and ApT = 1 GeV (for the other two ones). The J/ip decay branching ratio 
into the e~^e~ mode is also included. The magnitude of asymmetries and expected er- 
rors allows one to distinguish between different parametrizations of polarized parton 
distribution functions. Fig.4b shows expected asymmetries depending on the pseu- 
dorapidity of J/ip 8.t pr = 2 GeV. Statistical error bars correspond to the integrated 
cross section over bins At] = 1 and APt = 0.5 GeV. As in the previous case, the 
errors are small in a wide range of pseudorapidity interval and give a possibility to 
distinguish between different parametrizations for Ag{x). It is worth mentioning that 
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4 4.5 

PT(GeV) 

Figure 4. The expected asymmetries at ^/s = 39 GeV. Solid line corresponds Gehrmann- 
Stirling polarized parton parametrization, set A; dashed line, set C [|^. (a) expected 
asymmetries versus transverce momentum, (b) versus pseudorapidity of J/tp- Dot-dashed 
curve in the fig. a corresponds to the color-singlet contribution to asymmetry. 



gluon-quark collision subprocesses give about 5 -=-10% in Act. The dot-dashed curve in 
fig.4a represents the color-singlet contribution to the asymmetry for the parametriza- 
tion set A (only the color singlet contribution in Aa was taken into account). It is 
clear that the color octet contribution dominates in the expected asymmetries of J/i/j 
production. One may ask, where the such a difference in the shapes of asymmetries 
is coming from, as the subprocess asymmetries are not sensitive to colour structure. 
The main reason for this is the large matrix element < 0|08^'^(^S'o)|0 >. The pro- 
duction of such a state has a large asymmetry ^4 ~ 1 (at ~ 1, which is not far from 
the reality for the considered pr), providing the smooth transition to 2 — 1 process, 
discussed above. As such a state is absent in CSM, the large asymmetry itself is 
already a sign of a presence of colour-octet contribution. 

One of main parameters of the model is the mass of the charm quark. Fig. 5 shows 
the asymmetries of J/ ip production dependening on ttlc for two values of pt for a large 
set of polarized parton parametrizations (set A). As displayed in fig.5, the expected 
asymmetries are practically insensitive to the quark mass above m-c = 1.5 GeV. 
Therefore, the double spin asymmetry of J /if) production, unlike the cross section, 
should be free form uncertainties caused by a unknown mass of the intermediate color 
octet states. 

For comparison we have also calculated the expected double-spin asymmetries of 
J/ip production at RHIC energies. The results are given in Fig. 6 for two different 
values of energy. The expected statistical errors are calculated with the anticipated 
the integrated luminosities at the corresponding energies. As we can see from Fig. 6, 
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Figure 5. Expected double-spin asymmetries at y/s = 39 GeV versus mass of charm quark. 



the expected asymmetry decreases with increasing c.m.s. energy. The statistical 
errors are calculated by integration over pt with bins ApT = 0.5 GeV of the differential 
cross sections. 



5 Conclusions 

In this paper we investigated the expected double spin asymmetries of heavy quarko- 
nium hadropro duct ion at HERA-A^. To deal with experimentally observed quantities, 
we considered J/ip meson production at nonzero transverse momenta, pt > 1-5 GeV. 
Unlike the calculations of ||I6|, where only the lowest order subprocesses were taken 
into account (2^1), we considered J/ip production in the subprocesses 2 — > 2 be- 
cause large values of pt can not be caused by internal motion of partons. We have 
calculated the heavy quark pair color octet and color singlet S and P states produc- 
tion cross sections for different helicities of colliding partons. The FORTRAN codes 
for calculated cross sections are available by E-mail. For calculation of the expected 
hadronic asymmetries we used more reliable values for the color-octet long distance 
matrix elements, which are in good agreement with those extracted from the J/ip 
electroproduction data . 

The magnitude of expected asymmetries and statistical errors at HERA-A^ allows 
one to distinguish between different parametrizations for polarized parton distribution 
functions (Gehrmann and Stirling, set A and C). On the other hand, measuring the 
asymmetry would give a possibility to extract information about the color-octet long 
distance matrix elements and to check universality of factorization. We also calculated 
the J/ ip production asymmetries at RHIC energies. By comparing the magnitudes of 
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Figure 6. Expected asymmetries and statistical errors at the RHIC for two different ener- 
gies. 

the expected asymmetries at HERA-A^ and STAR, it becomes clear that the energy of 
the fixed target experiment is more preferable for the investigation of the charmonium 
production asymmetry. 
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